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Abstract 

We present a new derivation of the spectral gap of the totally 
asymmetric exclusion process on a half-filled ring of size L by 
using the Bethe Ansatz. We show that, in the large L limit, the 
Bethe equations reduce to a simple transcendental equation in- 
volving the polylogarithm, a classical special function. By solving 
that equation, the gap and the dynamical exponent are readily 
obtained. Our method can be extended to a system with an ar- 
bitrary density of particles. 
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1 Introduction 

The asymmetric simple exclusion process (ASEP) is a model of driven dif- 
fusive particles on a lattice with hard-core exclusion (for general review see 
Spohn 1991). The ASEP appears as a minimal building block in a large 
variety of models for hopping conductivity (Richards 1977), polymer rep- 
tation (Widom et al. 1991), traffic flow (Schreckenberg and Wolf 1998) or 
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surface growth (Krug 1997). In particular, the ASEP in one dimension is a 
discrete version of the Kardar-Parisi-Zhang (KPZ) equation (Halpin-Healy 
and Zhang 1995). In biophysics, the ASEP has been used to describe the 
diffusion of macromolccules through narrow vessels (Levitt 1973) and the 
kinetics of protein synthesis on RNA (MacDonald and Gibbs 1969); more 
recently, a mapping between sequence-alignment and the exclusion process 
has been proposed (Bundschuh 2002). Prom a theoretical point of view, the 
ASEP plays the role of a paradigm in non-equihbrium statistical mechanics: 
it displays a variety of features such as boundary induced phase transitions 
(Krug 1991), spontaneous symmetry breaking in one dimension (Evans et al. 
1995) and dynamical phase separation (Evans et al. 1998). 

Exact results for the ASEP in one dimension have been derived using 
two complementary approaches (for a review see Derrida 1998). The Matrix 
Ansatz (Derrida et al. 1993) allows to calculate steady state properties such 
as invariant measures (Speer 1993), current fluctuations in the stationary 
state and large deviation functionals (Derrida et al. 2003). The Bethe Ansatz 
(Dhar 1987) provides spectral information about the evolution (Markov) op- 
erator (Gwa and Spohn 1992; Schiitz 1993; Kim 1995) which can be used 
to derive large deviation functions (Derrida and Lebowitz 1998; Derrida and 
Appert 1999; Derrida and Evans 1999). The exact relation between these two 
techniques is still a matter of investigation (Alcaraz et al. 1994; Stinchcombe 
and Schiitz 1995; Alcaraz and Lazo 2003). 

The relaxation time to the stationary state for the ASEP on a lattice of 
size L scales typically as a power law, L^, z being the ASEP dynamical ex- 
ponent. The calculation of z for a one dimensional system by Bethe Ansatz 
is an important exact result that was first announced by Dhar, who found 
2; = 3/2 (Dhar 1987). The spectral gap (and thus z) was subsequently calcu- 
lated for the half-filling case by Gwa and Spohn (1992) and for an arbitrary 
density by Kim (1995) who mapped the ASEP into a non-Hermitian XXZ 
Heisenberg spin chain. The one-dimensional ASEP belongs to the KPZ uni- 
versality class and therefore the KPZ dynamical exponent in one dimension 
is equal to 3/2; this result was previously deduced from Galilean invariance 
and renormalization group arguments (for a review see Krug 1997). 

In this work, we present a new method of calculating the spectral gap of 
the totally asymmetric exclusion process (TASEP) starting from the Bethe 
Ansatz equations. Our method, based on an analytic continuation formula, 
circumvents the technical difficulties involved in the derivation of Gwa and 
Spohn (1992) and renders the calculation of the ASEP dynamical exponent 
much more concise and transparent. Besides, our technique can be readily 
extended to the arbitrary density case and allows us to calculate the spectral 
gap for the asymmetric exclusion process with a tagged particle. 
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The outline of this paper is as follows. In section 2, we recall the defi- 
nition of the TASEP, present the Bethe Ansatz equations without deriving 
them and summarize their analysis. In section 3, we present our original 
calculation of the spectral gap in the half-filling case. Concluding remarks 
and generalizations of our method are given in section 4. 



2 Bethe Ansatz for the TASEP 

2.1 The TASEP model 

We consider the totally asymmetric simple exclusion process on a periodic 
one dimensional lattice with L sites (sites i and L + i are identical). In 
this model, the total number n of particles is conserved. Each lattice site i 
(1 < i < L) is either occupied by one particle or is empty {exclusion rule). 
Stochastic dynamical rules govern the evolution of the system: a particle on 
a site i at time t jumps, in the interval between times t and t + dt, with 
probability dt to the neighbouring site i -|- 1 if this site is empty. The total 
number of configurations for n particles on a ring with L sites is given by 
Q = L!/[n!(L — n)!]. In the stationary state, all configurations have the same 
probability 1/n (Derrida 1998). 

A configuration can be characterized by the positions of the n particles 
on the ring, {xi,X2, ■ ■ ■ ,Xn) with 1 < xi < X2 < ■ ■ ■ < Xn < L. We call 
ipt{xi, . . . , Xn) the probability of this configuration at time t; the probability 
distribution ipf of the system at time t is thus a Q-dimensional vector. As 
the ASEP is a continuous-time Markov {i.e., memoryless) process, the time 
evolution of ipt is determined by the master equation 

where the transition rate Q x Q matrix M is the Markov matrix. A right 
eigenvector t/j is associated with the eigenvalue of M if 

Mi} = Eij. (2) 

The Markov matrix M is a real non-symmetric matrix and, therefore, its 
eigenvalues (and eigenvectors) are either real numbers or complex conjugate 
pairs. The spectrum of M contains the eigenvalue _E = and the associated 
right eigenvector is the stationary state iIj{xi, . . . ,a;„) = 1/^2. Because the 
dynamics is ergodic {i.e., M is an irreducible and aperiodic Markov matrix), 
the Perron-Probenius theorem implies that all eigenvalues E except have 
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a strictly negative real part; the relaxation time of the corresponding eigen- 
mode is r = — 1/Re(i?). (The imaginary part of E gives rise to an oscillatory 
behaviour) . 

In this paper, we shall calculate the gap Ei, i.e., the non-zero eigenvalue 
of M with largest real part. The eigenmode associated with Ei has thus the 
longest relaxation time that scales as L^, z being the dynamical exponent. 

2.2 The Bethe equations 

Writing M explicitly, the eigenvalue equation Q becomes 

i 

(3) 

where the sum runs over the indexes i such that < Xi — 1, i.e., such 
that the corresponding jump is allowed. The Bethe Ansatz assumes that the 
eigenvectors ip can be written in the form 

^l;{x,, . . . , x„) = ^ z^^l^ zll^-^ . . . zll^^ , (4) 

where E„ is the group of the n! permutations of n indexes. The coefficients 
{Ao-} and the fugacities {^i, . . . , z„} are complex numbers to be determined. 
The eigenvalue E associated with an eigenvector of the form (jH) is given by 

n 

E = -n + Y. ^I^i ■ (5) 

i=l 

Using matching conditions at the boundary surfaces = Xj — 1 and the 
periodicity of the lattice, it can be shown that a vector if) of the type (jH) is 
an eigenvector of M if the fugacities {zi, . . . , Zn} satisfy the Bethe equations 

n 

(^,-ir^-^ = -n(l-^.) for 1 = 1,. ..,n. (6) 
i=i 

The procedure for deriving these equations has been thoroughly explained in 
Halpin-Healy and Zhang (1995) and Derrida (1998). 

The obvious solution 2:1 = . . . = z„ = 1 provides the stationary distri- 
bution with eigenvalue 0. More generally, given a solution {zi,...,Zn\ of 
Eq. the corresponding eigenvalue E is obtained from Eq. (0); moreover, 
the coefficients {A„} and thus the eigenvector ip are uniquely determined. 
In order to have a complete basis of eigenvectors, VL independent solutions of 
the Bethe equations (jHl) are needed. 
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Following Gwa and Spohn (1992), we introduce Zi = 2/zi — 1. The 
equations (0) and ^ then become, respectively, 

2E = -n + J2Z,, (7) 

and 

(1 - (1 + Z,)^-" = -2^ n 1^ fo^ z = l,...,n. (8) 

We notice that the left hand side of Eq. (jHl) is a polynomial in Zi whereas 
the right hand side (r.h.s.) is independent of the index i. 

The analysis of the Bethe equations is simplified if only half-filled models 
are considered, that is, if L = 2n (Gwa and Spohn 1992). Equation (jH} then 
reduces to 

n- 2 — 1 

(l-Z^r = -rl[^— for z = l,...,n. (9) 

The half- filling restriction does not affect the physical behaviour of the ASEP: 
in the large L limit, models with arbitrary density, p = n/L G ]0, 1[, belong to 
the same universality class. However, systems with vanishingly small density 
of particles [p — *■ 0) or holes (p — > 1) exhibit a different behaviour and will 
not be discussed here. 

2.3 Analysis of Bethe equations 

Taking advantage of the fact that the r.h.s. of Eq. is independent of the 
index i, the Bethe equations can be reformulated as follows. Consider the 
polynomial equation 

(i-z2)" = y, (10) 

where F is a given complex number. Writing 

F = -e"", (11) 
u being a complex number with —1 < Im(-u) < 1, we obtain the n-th roots 

of y 

= Q{u+i)n/n^{m-l)2in/n for m = 1 , . . . , n . (12) 

The ?/m's are evenly spaced on a circle of center and radius and are 

labeled counter-clockwise < Arg(?/i) < Arg(?/2) < . . . < Arg(|/„) < 27r. 
Thus, the 2n solutions {Zi, . . . , Z2n) of Eq. (fTUj) are 

Zm= {I- ymf^ ; Zm+n = -Zm with m = l,...,n. (13) 



O. Golinelli, K. Mallick — ASEP spectral gap 



6 



The branch cut of the function z^^^ is, as usual, the real semi-axis (— oo,0], 
i.e., for m = l,...,n, the argument of belongs to [— 7r/2, 7r/2[. We 
explain in Appendix |Bl that each is an analytic function of Y in the 
complex plane with a branch cut along [0, +00) and that the locus of the 
Zm's is a remarkable curve called a Cassini oval. 

We now choose n different roots (Zc(j))j=i^n among (Zi, . . . , Z2„), such 
that the choice function c:{l,...,n} —*>{!,..., 2n} satisfies 

1 < c(l) < . . . < c(n) < 2n. (14) 

There are precisely Q such choice functions, Q = (2n)!/n!^ being the size of 
the Markov matrix. Finally, we define 

^(F) = -4" n , (15) 

where is a function of Y and of the choice function c. The Bethe equa- 
tions (jni) are equivalent to the self- consistency equation 

Ac{Y) = Y. (16) 

Given the choice function c and a root Y of this equation, the -^c(j)'s are 
determined from Eq. ()10j) and the corresponding eigenvalue is obtained 
from Eq. ((Tj). 

For small values of n, the above described procedure allows us to compute 
numerical solutions of the Bethe equations. From our numerical observations, 
we conjecture that for each choice function c (among the Q possible choice 
functions), the self-consistency Eq. (fTBj) has a unique solution Y that yields 
one eigenvector ip^ and one eigenvalue E^. This suggests that the Bethe 
equations yield a complete basis of eigenvectors for the ASEP. 

Let us first consider the choice function c(j) = j, i.e., the selected Z/s 
are the n solutions of Eq. (|Tn|l with largest real parts. As this choice plays 
an important role in the following analysis, we define 

MY) = -^-"fi^^, (IT) 

n 

2Eo = -n + Y,Zj. (18) 
i=i 

We emphasize that Eq is an implicit function of Y. The equation ^o(^) = Y 
has the solution Y = that yields Zj = 1 for all j and provides the stationary 
distribution (or ground state) with eigenvalue 0. 
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Numerical observations (Gwa and Spohn 1992) indicate that the first 
excited eigenvalue Ei corresponds to the choice c(j) = j for j = 1, . . . , n — 1 
and c(n) = n + 1; i.e., the first excited state is obtained from the ground 
state by the excitation (n — n + 1) . Writing Ai{Y) and Ei for the functions 
Ac and Ec corresponding to this choice function, we have, from Eqs. (fTHj fTTI 
andHH), 

2E, = 2Eo-(Zi + Z„), (20) 

where we have used ^„+i = —Zi. The excitation (1 2n), that is, c(j) = 
j + 1 for j = 1, . . . , n — 1 and c(n) = 2n, also leads to Eqs. (fT^ and 
and thus to the same eigenvalue Ei. The first excited state has therefore a 
degeneracy of order 2. 

Consequently, in order to find the expression for the gap Ei, we must 
solve the self-consistency equation 

AiiY) = Y, (21) 

calculate the Z^-'s for j = 1, . . . , n, and finally deduce Ei from Eq. (pUj) . 

In the above discussion, we closely followed Gwa and Spohn (1992) to 
present the Bethe Ansatz equations for the TASEP. We shall now solve these 
equations and calculate the gap by a radically different and simpler method. 



3 Calculation of the gap 

Let us define E{Y) as 

Ao(r) = rexp(F(y)). (22) 

In Appendix we have derived the following identities, valid for \Y\ < 1, 

FiY) = t^y, (23) 
k=i 

oo 

-AE, = Y^'^Y', (24) 
fc=i 



where the w^'s are given by 
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From the Stirling formula, the leading order behaviour of Wk for A; ^ oo, is 
given by 

Wk-^^. (26) 

From the power series and we deduce that ^o(^) and Eq are analytic 
functions of the complex variable Y inside the unit circle. This property 
is not obvious a priori: the functions Aq(Y) and Eq, defined in Eqs. (fT7|) 
and (fTHj) . respectively, depend implicitly on Y via the Z^s that involve a 
branch cut along [0, +oo). Indeed, for a generic choice function c(j), Ac(Y) 
and Ec are analytic only in the complex Y plane with a cut along [0, +oo) 
and, therefore, are not analytic in the neighbourhood of F = 0. This special 
property of ^o(^) and Eq is obtained by an explicit calculation in Appendix 
1X1 and from geometrical considerations in Appendix ^ 

Using Eqs. (fT^ . and (j^ . the self-consistency equation (pij) . that 
determines the gap, reduces to 



oo 



From Eq. and Eq. ((211), we have 

Wkn-l ^^k 

k=l 

Combining Eq. (fTfj) and Eq. (j^ . we obtain 



In T-/ . (27) 



4^1 = -^r^ + 2^1 + 2^„ . (28) 



\k{2kn-l) V l + ^i/ V 1 + Z 



(29) 

where we have used {2kn)wkn = {^kn — l)wkn-i- 

Thus, to find the gap for a half-filled system with n particles, we must 
solve Eq. ^I7\i for Y and substitute the result in Eq. (|^. We emphasize that 
the power series in these equations represent, inside the unit disk, analytic 
functions of Y that are defined in the whole complex plane with a cut along 
[l,+oo). 

We now consider the thermodynamic limit, n — >■ oo. We obtain, at leading 
order, from Eq. 



oo 



FiY) = Y.^Y' - -^^-^s/2in (30) 



k=l 



OO "1 
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where we have used the polylogarithm function of index s 



LiJz) 



oo 



E 

fc=l 



(32) 



By virtue of the integral representation 



LiJz) 



T(s) 



dt 



(33) 



the function Li^ can be extended by analytic continuation to the whole com- 
plex plane with a branch cut along the real semi-axis [1, +cxo). In the large 
n limit, we deduce from Eqs. pil IT^ andlT ^ that 



n 



n 



3/2 



(34) 
(35) 



where we have supposed that Y = —e^^ remains finite when ri — > oo. Using 

i-z _ 

l+Z ^ 3' 



these expressions and the expansion In j— ^ = — 2Z — + 0{Z^), Eq. (j2 



reduces to 



Li3/2(-e"") = -271 



-u + i 



,1/2 



—u — I 



and the gap (pUj) . at the leading order, is given by 



El 



n- 



3/2 



Li 



5/2 ( 



TT 



3/2 



6 



{-u + if'^ + {-u - i) 



N3/2 



(36) 



(37) 



[Notice that the r.h.s. of Eq. ()36|) and of Eq. ()37|) are real when u is real.] 
With the help of the Maple software, we find a unique solution of Eq. fj36|l 
in the strip — 1 < Im(ti) < 1 that is real and is given by 



u = 1.119 068 802 804474. . . 



(3^ 



Inserting this value of u in Eq. ()37|1 yields the large n (or large L) behaviour 
of the gap 



El 



2.301345 960 455 050... 



n 



3/2 



6.509 189 337976 136... 

Z372 ■ 



(39) 



This is precisely the result obtained by Gwa and Spohn (1992). This gap 
scales as L~^/^ and is real for the TASEP in the half-filling case. 
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4 Summary and discussion 

In this work, we have calculated the gap of the TASEP in the limit of a 
large size system by using the Bethe Ansatz. We first take the large n 
limit of the Bethe equations inside the unit circle, then perform the analytic 
continuation of these equations in the whole complex plane with a branch 
cut along [l,+cx3) and finally solve them. Gwa and Spohn (1992), on the 
contrary, first represent the analytic continuation of the Bethe equations for 
a fixed value of n as a n-dependent complex integral (thanks to the Euler- 
Maclaurin formula) and then extract the gap from the large n limit of this 
integral representation which is rather a delicate operation. We have shown 
here that the derivation of the TASEP gap is greatly simplified by performing 
the calculations in the reverse order, that is taking the large n limit first and 
the analytic continuation afterwards. 

We do not claim that it is always true that large n limit and analytic 
continuation are commuting operations. If the solution Y of Eq. (j?T|) . in the 
large n limit, diverges to oo or approaches asymptotically the branch cut, 
reversing the order of operations may not be possible. Fortunately, in our 
problem, Y remains a bounded negative real number when n — oo. 

We have applied our method to several other problems but, for sake of 
conciseness, we simply list these additional calculations without giving the 
details. Our technique provides the subleading corrections to the gap and 
allows us to calculate the eigenvalue of the highest excited state, of a finite 
excitation above the ground state, or below the highest excited state. The 
gap for the totally asymmetric exclusion process with an arbitrary density p 
can also be obtained by elementary means, and we find in agreement with 
(Kim 1995) 

E, (p) = 2^p{l - p) E, (p = 1/2) ± _(2p - 1) . (40) 

where Ei{p = 1/2) is the gap for the half-filling case given in Eq. p9|l : we 
notice that for a density other than one-half the gap has a non-zero imaginary 
part. We have also studied generalizations of the ASEP by introducing a 
tagged particle that has the same dynamics as the other particles: the gap 
then scales as 

Finally, we emphasize that the formula ()33|1 already appeared in the work 
of Derrida and Appert (1999): indeed, Eqs. (fTfjl and are similar to those 
used in their calculation of large deviation functions of the ASEP by Bethe 
Ansatz. 
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A Derivation of Eq. (ESI) and Eq. (El) 

The numbers Wk defined in Eq. fl25|l are the coefficients of the Taylor series 

-1 oo 

-^— = T.^kx'. (41) 



By integration, we find 



1 ^y^k-i. 



VT^ = i-^E^xK (42) 

Recalling that Zm = v^l — ym [see Eqs. (|T^ and |T!?jl ]. we obtain 

n 1 oo n 

T.Zm = n-l y^T^yt- (43) 

m=l ^ k=l m=l 

The fact that the ?/m's are the n-th roots of Y leads to the following relation 

is a multiple of n , ^ 

^"^ ~ 1 otherwise. ^ 

r7i=l 

Inserting this relation in Eq. (j43|) . we obtain 

" 1 oo 

m=l ^ fc=l 

We thus find, thanks to the crucial identity (jH)), that Y^^^i is analytic in 
Y inside the unit circle. Finally, substituting Eq. (j45|l in Eq. (fT8|l . we obtain 
Eq. (121. 

The derivation of Eq. follows similar steps. We first notice that the 
Taylor expansion of the function 



is given by 



= (46) 



m = T.T^'- (47) 

k=l 
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[This follows from /(O) = and f'{x) = \ (^-^ - l) ]. Using Eq. ^ and 
the identity 11^=1 = deduce that 



n 



4 I- Z„ 



From Eq. ^ and Eqs. (HHI), and (gHl), we obtain 



m=l k=l m=l 

This equation reduces to Eq. by virtue of the identity 



B Roots of the Bethe Equations and Cassini 
ovals 

The polynomial equation 

(1 - Z^)" = Y (50) 

where F is a fixed complex number, has 2n solutions, {Zi, . . . , Z2n)- The 
purpose of this appendix is to explain how these solutions can be labeled in 
a coherent way so that each root ZmiX) is an analytic function of Y . 

We first notice that defined in Eq. flT^ . is an analytic function of 
Y in the complex plane with a branch cut along the real semi-axis [0, +cx3). 
Nevertheless when Y crosses [0, -|-cx)), the functions yi,y2 ■ ■ ■ ,yn are the ana- 
lytic continuations (above the axis) of respectively yn,yi, ■ ■ ■ ,yn-i- thus the 
existence of the branch cut along [0, +oo) is due to the labeling of the roots. 

The complex numbers (Zi, . . . , Z2n) belong to the curve defined by 

\Z-l\\Z + l\=r with r = |F|i/" (51) 

and called a Cassini oval. A Cassini oval is the conformal transformation 
of the circle of center 1 and radius r by the function z — ^ z^^'^. Its shape 
depends on whether the point is inside or outside the circle, i.e., whether 
r < 1 or not. When r < 1, the curve of equation (jKT|) consists of two ovals 
around the points Z = ±1. The numbers [Zi, . . . , Z„) lie on the right oval 
and {Zn+i, . . . , Z2n) on the left oval. For the marginal case, r = 1, the curve 
is the lemniscate of Bernoulli, with a multiple point at Z = 0. When r > 1, 
the Cassini oval is a single loop with a peanut shape (when r g]1,2[) or an 
oval shape (r > 2). See Fig. Q where the cases r = 0.9, 1 and l.l are drawn. 
In the large-r limit, the oval tends to the circle of radius ^/r. 
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z, 



z, 



n-1 




n 



Figure 1: Labeling the roots Zm of the equation (1 — Z^)" = Y . Here Y = 
^i(t>^n yjj^if^ ^ — ]^g^ — 77/2, and r G {0.9, 1, 1.1}. The continuous curves are 
the corresponding Cassini ovals (see text). When r is fixed and (j) goes from 
to 2tt, each slips counter-clockwise along the Cassini ovals. Then, the 
jump (f) = 27T ^ 0, i.e., Y crosses [0,+oo), consists of a global shift of the 
labels m around each continuous curve. 

We now discuss the analyticity properties of ZmiY): Z^ is an analytic 
function of ym with a branch cut along (—00, 0]; this branch cut is compatible 
with that of ym{Y). Consequently Zm is an analytic function of Y with a 
branch cut along [0, +0x3). Moreover when Y crosses the real segment [0, 1[, 
the functions Zi{Y), Z2{Y) . . . , Z^iY) are the analytic continuations (above 
the axis) of the functions Zi{Y), . . . , respectively. See Fig.[TJ 

(A similar property is true for the functions Zn+i(Y), Z„+2(F) . . . , Z2n(Y) 
that belong to the left oval). But when Y crosses ]l,+oo), the functions 
Zi{Y), Z2{Y), . . . , Z2n(Y) are the analytic continuations (above the axis) of 
the functions Z2n{Y), Zi{Y), . . . , Z2n~i{Y) respectively. Consequently the 
branch cut of the function Aq{Y) defined in Eq. (fT7|) is [1, +00) and not 
[0,+oo). 
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